home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / CORREL.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  619b  |  23 lines
  1. PROCEDURE correl(data1,data2: glnarray; n: integer; VAR ans: gl2narray);
  2. (* Programs using routine CORREL must define the type
  3. TYPE
  4.    glnarray = ARRAY [1..n] OF real;
  5.    gl2narray = ARRAY [1..2*n] OF real;
  6. in the main routine. *)
  7. VAR
  8.    no2,i,ii: integer;
  9.    dum: real;
  10.    fft: gl2narray;
  11. BEGIN
  12.    twofft(data1,data2,fft,ans,n);
  13.    no2 := n DIV 2;
  14.    FOR i := 1 TO (no2+1) DO BEGIN
  15.       ii := 2*i;
  16.       dum := ans[ii-1];
  17.       ans[ii-1] := (fft[ii-1]*ans[ii-1]+fft[ii]*ans[ii])/no2;
  18.       ans[ii] := (fft[ii]*dum-fft[ii-1]*ans[ii])/no2
  19.    END;
  20.    ans[2] := ans[n+1];
  21.    realft(ans,no2,-1)
  22. END;
  23.